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The review of geochronological and historical data documents that the largest southern European deltas 
formed almost synchronously during two short intervals of enhanced anthropic pressure on landscapes, 
respectively during the Roman Empire and the Littie Ice Age. These growth phases, that occurred under 
contrasting climatic regimes, were both followed by generalized delta retreat, driven by two markedly 
different reasons: after the Romans, the fall of the population and new afforestation let soil erosion in river 
catchments return to natural background levels; since the industrial revolution, instead, flow regulation 
through river dams overkill a still increasing sediment production in catchment basins. In this second case, 
furthermore, the effect of a reduced sediment flux to the coasts is amplified by the sinking of modern deltas, 
due to land subsidence and sea level rise, that hampers delta outbuilding and increases the vulnerability of 
coastal zone to marine erosion and flooding. 



Deltas are strongly three-dimensional lithosomes that represent one of the most complex, fragile and 
densely populated landscapes of mid-latitude modern coastal systems. Understanding the life cycle of 
modern deltas since their embryonic stages is crucial for disentangling the role of natural vs. anthropic 
forcing, an open issue that bears substantial implications on our full awareness on the management of the entire 
coastal zone. By underpinning the timing of the main phases of delta growth it is possible to define how 
pervasively have human populations modified landscapes and to predict threats to delta stability during forth- 
coming, and possibly abrupt, environmental changes. The growth of deltas reflects the forcing by eustatic 
oscillations, which determines the available space for sediment accumulation, the riverine input, defined by water 
and sediment discharge, and the oceanographic regime, responsible for sediment dispersal 1 . The lithology and 
morphology of the drainage basin and the climatic regime govern sediment production and delivery to the coast, 
while the power of marine processes (storms, waves and tides) and local vertical movements of the delta plain 
controls sediment accumulation and the final shape of the delta 2 . 

The evolution of modern deltas can be reconstructed combining complementary information from multi- 
proxy analyses, at core scale, and from geomorphological reconstructions of the growth of delta plains, channel 
avulsion and delta lobe progradations 3 . The comparative study of the Ebro, Rhone, Po and Danube deltas, 
representing the four main deltas of the northern Mediterranean and Black Sea regions (Fig. 1), is favored by 
the following circumstances: 1 - the oceanographic regime of the receiving basins has allowed delta formation and 
substantial preservation; 2- the four catchments present comparable climate conditions having contiguous 
watersheds (comprised between 40° and 48° N) impacted by human diffusion since the Neolithic; 3- the ages 
of the main progradation phases are constrained by robust chronologies, also supported by map documentations 
over the last 5 centuries, providing the possibility to establish reliable trans-delta time lines. The comparison of the 
evolution of northern Mediterranean and Black Sea deltas, permits a super-regional correlation of their growth/ 
retreat phases; this approach to the study of modern deltas allows understanding the relative role of climate 
forcing and of growing human population in their evolution. 

To identify the intervals of maximum rate of delta outbuilding and to attempt a chronologically constrained 
comparison of phases of accelerated delta growth among different systems, it is useful to simplify the three- 
dimensional progradation of each individual lobe and to convert it into "spatially-averaged" linear values 
(Supplementary SI), quantifying also the evolution of the paleoenvironmental constraints (climate, relative sea 
level and oceanography, Supplementary S2). In this way it is possible to identify synchronous changes in the rates 
of delta growth that may reflect anthropic perturbations of the natural landscapes, likely in response to the spread 
of new techniques in forest clearance, cropland cultivation, or urbanization. By comparing the linear rates of delta 
front progradation (Fig. 2, Supplementary S3) we differentiate the impact of Mid- Late Holocene high-frequency 
climatic oscillations from that of increasing human population on sediment production and, ultimately, delta 
construction and preservation. 
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Figure 1 | Location of the four delta systems with catchments and growth phases. Drainage basins of the largest southern European river 
systems on SRTM 90 m Digital Elevation Data (CGIAR-CSI), with major dams and reservoirs (red dots) and location of the Soppensee Lake 
(black star; see pollen spectra in Figure 3). The inset graph shows the number of river dams built per year since 1900 AD; the temporary decline in dam 
construction in the '40 can be attributed to the Second World War (see Supplementary S5). Below: satellite views of modern deltas (from the 
Landsat Archive and the Global Land Survey), with red lines representing their embryonic stages, during the last phases of post-glacial sea level rise 
(between 9 and 5.5 kyr BP), and the white lines representing the main dated progradational events of delta growth (see Supplementary S3). 



Results 

Delta growth, climate and humans. The onset of the northern 
Mediterranean and Black Sea deltas follows an ancestral phase 
dominated by estuary fill and relatively slow delta growth, that can 
be traced back to the Mid-Holocene (Fig. 2), around 6 kyr BP, when 
deceleration in the rate of sea level rise allowed sediment accu- 
mulation at the river mouths 4,5 . The following few thousand years 
of delta growth reflect a complex picture of high frequency climatic 
changes 6 accompanied by the spreading of human populations 7 , 
impacting the land cover, the hydrological cycle and ultimately the 
sediment production and transport through river catchments 
(Fig. 3). During the last 3 kyr, in particular, a strong relation was 
observed between the fall of past civilizations and environmental 
changes driven by climatic deterioration 8 . In Europe the first impor- 
tant human diffusion dates back to the Late Iron Age and continues 
until the first few centuries AD, as testified by the increased 
deforestation for construction activities 9 . Extensive tree ring-based 
reconstructions over central Europe show that this interval was cha- 
racterized by mild climatic conditions, with summer temperatures 
and precipitations substantially above average 9 . This evidence is 
supported by high values in Greenland surface temperature and by 



a persistent positive North Atlantic Oscillation (NAO) index 1011 . The 
rise of classical civilizations and the spread of agricultural activities 
led to an increase in population and to significant deforestation in 
central and western Europe 12 . The peak of human growth was 
reached at the apex of the Roman Empire around 250 yr AD 12 , 
which corresponds also to the first hemispheric-scale pollution 
event caused by metal extraction and production 13 . This important 
anthropic control on the natural environment at European scale 
resulted in a substantial increase of soil erosion and in the first 
major burst in the outbuilding of river- dominated deltas that 
interrupted, around 0 AD, the previous phase of slower strand 
plain construction. The expansion of human populations across 
central Europe during the six centuries of the Roman Empire was 
accompanied by marked modifications of landscapes, with forest 
clearance and land-use changes for cropland and pasture; these 
changes resulted in a huge increase in sediment production in the 
river watersheds that led to the first synchronous accelerated phase of 
Pan-Mediterranean delta construction (Fig. 2). During this interval 
several roman harbors were rapidly mud-filled in response to 
enhanced river run off and increased sediment load 14 . 

The main environmental signature of the fall of the Roman 
Empire, at the end of the 4 th century, is the new afforestation, also 
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Figure 2 | Correlation of linear progradation rates of the four delta systems. Summary of the main intervals of delta growth of the Ebro, Rhone, Po and 
Danube systems, reporting their maximum rates as orders of magnitude (see Supplementary S 1 and S3). Note that the time window is in two scales with 
1 kyr steps, before 1500 AD, and 20 yr steps, onwards; delta growth rates are also summarized against a linear time scale in Figure 3. Two major and 
substantially synchronous phases of delta construction correspond to the expansion from 300 yr BC to 250 yr AD, and the last ca. four centuries, 
matching the Little Ice Age. Retreat of major delta fronts occurred after the Roman Empire and impacted almost all delta lobes during the last 50 yr 
(red bars). 



simulated by the HYDE model for central Europe 12,15 ; this afforesta- 
tion reflects a substantial reduction in the number of inhabitants in 
response of the barbarian invasions during the Migration Period 
(5 th -9 th centuries AD), accompanied by an overall decrease of air 
temperatures, as indicated also by pollen spectra from the 
Switzerland Alps (Fig. 3). All the southern European deltas show a 
generalized retreat during the millennium following the decline of 
the Roman Empire; this interval encompassed repeated climate oscil- 
lations, including the cold Dark Ages and the Medieval Climate 
Anomaly (MCA) that was accompanied by a gradual socioeconomic 
consolidation with a consistent demographic growth. The MCA, a 
250-yr period of climatic amelioration, shows a closely-spaced suc- 
cession of warm and wet summers that favoured population growth 
and consequent forest clearance for cropland and pasture. NAO 
index reconstructions 16 show a prolonged positive state during the 
MCA 17 that favored a period of increased river discharge reflect- 
ing enhanced precipitations and retreat of Alpine glaciers 18 . One 



outcome of this trend was the major avulsion phase of Ficarolo in 
the Po River Delta ( 1 1 52 AD; Fig. 2 and 3 ) . 

A first cold spell around 1300 AD marked the transition from the 
MCA to the Little Ice Age (LIA, 1400-1850 AD) when the deteriora- 
tion of climatic conditions prepared the ground for the pandemic 
diffusion of the Black Death from (1347- 1350 AD), that caused a 60% 
reduction in European population that was accompanied by a new 
major phase of afforestation 12 . The European population returned to 
pre-pandemic conditions only around the second half of 1400 AD 
and has increasingly grown ever since. 

The LIA marked a dominant reorganization of large-scale atmo- 
spheric and oceanic circulation patterns, as highlighted by a typically 
negative and more variable NAO index 1116 . In central and southern 
Europe, climatic deterioration entailed a strong reduction in temper- 
ature and precipitation, with a decrease in winter temperature up to 
0.5°C, that led to several phases of Alpine glacier expansion and 20- 
50 mm negative anomalies in precipitation 1619 , with consequent 
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Figure 3 | Multi-proxy reconstructions of climate change, human diffusion and delta progradation. The last 3 millennia were characterized by several 
climatic oscillations, synchronous at northern Hemisphere scale, with two major extremes in the Medieval Climate Anomaly (MCA) and the Little Ice Age 
(LIA). The trend of pollen record from the central Alps 30 reflects temperature and precipitation changes all over the Europe region' and can be compared 
with the temperature record derived from Greenland ice core 10 . Growth rates of the main four northern Mediterranean and Black Sea river deltas (see 
Supplementary SI, each step is one order of magnitude) increased during the Roman Empire, and dramatically since the onset of the Little Ice Age (LIA), 
characterized by closely-spaced intervals of advancing and waning of Alpine Glaciers 19 , increased flood frequency 21 , high storm activity 31 , but only a slight 
increase of total precipitation 9 . These natural factors may reflect the occurrence of repeated and closely spaced NAO-index anomalies ( n , positive, red, 
and negative, blue) and several centuries of overall decreasing T 10 . While human expanded, in the past, during intervals of milder climate (Roman 
Empire), this was the first substantially cooler interval severely impacted by humans through, inter alia, forest clearance 12 , increasing population 12 and 
related land consumption. The change in land use is well reflected by increased shrub + herb vs. tree in pollen spectra both in the alpine region and in the 
Mediterranean, especially in the second half of the LIA. 
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Figure 4 | Modern water and sediment flux. Coupled diagram of suspend sediment load (brown) and fresh water (light blue) demonstrating the 
impact of anthropic reservoirs on river discharge in all four Mediterranean deltas; gray stripes highlight the timing of construction of major dams 
(Supplementary S3 and S5). Peaks of water discharge are in phase with peaks in suspended sediment load only before dam's construction until the '60. 




changes in soil vegetation cover 20 . By the second half of the 17 th 
century, roughly corresponding to the beginning of the Maunder 
Minimum (1645-1715 AD), a new period of delta growth started 
and continued for the entire LIA and 50 yr beyond. The climatic 
conditions during the Maunder Minimum, the coldest phase of the 
LIA, seem to be similar, although less amplified, to those during the 
Last Glacial Maximum 21 . Interestingly, the LIA had less dramatic 
impact on human societies with respect to previous intervals of cli- 
matic deterioration, because the population had become less vulner- 
able to climatic fluctuations thanks to technological developments, 
urbanization and industrialization. In this period, the growth of 
human population is the main driver for anthropogenic changes in 
land use and land cover at the expenses of natural forests. Hu- 
man population has grown exponentially since the XIX century, 



demanding more usable land for food, fuel and construction; all 
together, these factors lead to a generalized land-use change at 
European scale. The resulting regional modifications in hydrology 
and climate, through changes in albedo, surface roughness, runoff 
and groundwater storage 22 , determined a dramatic increase in sedi- 
ment production, sustaining a major and synchronous prograda- 
tional phase of all the southern European deltas (Fig. 2 and 3). This 
extremely enhanced sediment production due to anthropic activities 
was also reinforced by the occurrence of major flood events, con- 
nected to increased storm frequency and Alpine glacier collapse 23 . In 
the interval between 1650 and 1900 AD, in particular, sediment sup- 
plied by all rivers allowed substantial construction of delta plain in 
each system, and up to 350 km 2 for the Po River delta, with rates of 
linear basinward delta progradation up to 250 m/yr (Fig. 1 and 2). 
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The last century: deltas and the industrial era. Since the beginning 
of the Industrial Revolution, human control on the natural 
environment became overwhelming: at present one-third to one- 
half of land surface has been impacted and transformed by human 
activities 24 , with dramatic perturbations of global ecosystems. The 
anthropic impact on rivers and riparian zones, mainly for energy and 
irrigation, was achieved by fragmenting the river courses through 
dams and channel diversions 25 , that led to a dramatic decrease in the 
supply of coarse-grained material to the sea (Fig. 4) 26 . Due to 
sediment retention on anthropic reservoirs, the reduced riverine 
clastic input no longer balances the sink of delta plains caused by 
subsidence and sea level rise 27 . Delta systems become, therefore, 
more vulnerable to the action of marine processes, leading to a 
generalized phase of delta erosion and retreat. During the last five 
decades, several km 2 of delta plains have been submerged and all 
northern Mediterranean and Black Sea deltas underwent extensive 
erosion (Supplementary S4). The same fate has been documented 
for deltaic systems worldwide where similar human impact has 
occurred 28-29 . 

Discussion 

The quantitative comparison and correlation of the growth phases of 
the largest southern European deltas show that their outbuilding 
peaked almost simultaneously in two short intervals when human 
activities were at their maximum: the Roman Empire and Little Ice 
Age. These growth phases occurred during contrasting climatic con- 
ditions (dominantly warmer during the Roman Empire and cooler 
during the LI A). In both cases the anthropic control on the natural 
environments overshot the climatic signal forcing delta systems 
toward a man-made river-dominated signature. Both phases of 
man-made delta outbuilding were followed by discrete intervals 
of delta retreat caused by a dearth in sediment flux and, conse- 
quently, an increased effect of marine reworking. However, the two 
retreat phases, though similar in their impacts on deltas, have a 
markedly different origin: after the fall of the Roman Empire soil 
erosion was contrasted by a new afforestation phase and delta sys- 
tems just relaxed to their background (pre-anthropic) natural 
equilibrium; on the contrary, during the last 50 years, Medi- 
terranean and Black Sea deltas reflect the simultaneous increase in 
sediment production through soil erosion but yet reduced sediment 
flux to the coasts because of retention within reservoirs, fluvial regu- 
lation and excavation. The chronological constraints available for 
Mediterranean and Black Sea deltas indicate how natural systems 
respond "instantaneously" to human perturbation; this rapidity of 
response should be taken seriously into account in any management 
attempt of the coastal zone. 

Methods 

The data presented in this Manuscript are derived from a systematic scrutiny of up to 
50 scientific papers and reports on Mediterranean and Black Sea deltas, available in 
national and international literature (Supplementary S3). The modern southern 
European deltas were intensively investigated during the last two decades by integ- 
rating subsurface (seismic sections and boreholes) and subaerial data (satellite images 
of delta plains), supported by sedimentological, paleoenvironmental and chronolo- 
gical constraints. The main phases of delta progradation derived by the literature were 
constrained chronologically by hundreds of radiocarbon and few tens of OSL dates. 
The resulting chronologies are integrated for the last few centuries by historical 
records, in the form of ancient maps and chronicles, that allowed independent esti- 
mates of the main constructional and retreat phases in each delta. 

The averaged linear progradation rates presented here are quantified considering 
the three-dimensional nature of each delta system (Supplementary SI), and taking 
into account all the available regional paleoenvironmental constraints (climate, 
relative sea level oscillations, and oceanography, Supplementary S2). 
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